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Photoproduction reactions occur when the electromagnetic field of a relativistic heavy ion interacts 
with another heavy ion. The STAR collaboration presents a measurement of p and direct 7r + 7r~ 
photoproduction in ultra-peripheral relativistic heavy ion collisions at v /sjvjv=200 GeV. We observe 
both exclusive photoproduction and photoproduction accompanied by mutual Coulomb excitation. 
We find a coherent cross-section of cr(AuAu — > Au'Au* p°) = 530 ± 19 (stat.) ± 57 (syst.) mb, 
in accord with theoretical calculations based on a Glauber approach, but considerably below the 
predictions of a color dipole model. The p° transverse momentum spectrum ( p% ) is fit by a 
double exponential curve including both coherent and incoherent coupling to the target nucleus; we 
find (Tinc/o" C oh = 0.29 ± 0.03 (stat.) ± 0.08 (syst.). The ratio of direct tt + tt~ to p production is 
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comparable to that observed in collisions at HERA, and appears to be independent of photon 
energy. Finally, the measured p° spin helicity matrix elements agree within errors with the expected 
s-channel helicity conservation. 

PACS numbers: 25.20.Lj, 13.60.-r 



I. INTRODUCTION 

In heavy ion collisions when the electromagnetic field 
of one nucleus interacts with another nucleus, photopro- 
duction can occur P, Photoproduction is visible in 
ultra-peripheral collisions (UPCs), which occur when the 
impact parameter b is more than twice the nuclear radius 
Ra, so no hadronic interactions are present. The elec- 
tromagnetic field of a relativistic nucleus may be repre- 
sented as a flux of almost-real virtual photons, following 
the Weizsacker- Williams method [3j. In this framework, 
the physics of the interactions between particles is equiv- 
alent to that of the interactions between photons and 
particles. The photon flux scales as the square of the 
nuclear charge and so the cross-sections can be large in 
heavy ion interactions. 

Photoproduction of p°s occurs when a photon from one 
nucleus fluctuates to a quark-antiquark pair, which then 
scatters elastically from the other nucleus, emerging as a 
p° . The elastic scattering can be treated as being due to 
Pomeron exchange Q. 

The cross-section for p° production depends on the qq 
coupling to the nuclear target. For p° production at large 
transverse momentum, pt, the qq pair couples to the 
individual nucleons. The incoherent cross-section scales 
roughly as the atomic number A, minus a correction due 
to nuclear absorption of the p°. 

At smaller pt, roughly < h/ Ra-, the qq pair couples 
coherently to the entire nucleus; naively, this leads to a 
cross-section that scales as A 2 . The coherent produc- 
tion is regulated by the nuclear form factor F(t), so p° 
photoproduction is sensitive to the /9°-nucleon interaction 
cross-section and the nuclear structure functions 

There are three published calculations of the coherent 
p° photoproduction cross-section in heavy ion collisions. 
The first model (Klein and Nystrand - KN), uses vec- 
tor meson dominance (VMD) plus a classical mechanical 
Glauber approach for nuclear scattering. KN uses in- 
formation from the 7p — > Vp experiments for extrapola- 
tion Q . The model predicts a total coherent p° photopro- 
duction cross-section cr p o=590 mb in gold-gold collisions 
at y/sNN=200 GeV. The second model (Frankfurt, Strik- 
man and Zhalov - FSZ) treats the p° production using 
the generalized quantum VMD and the QCD Gribov- 
Glauber approach @, 0|; it predicts cr p o = 934 mb Q, 
about 60 % higher than the KN model, but with a sim- 
ilar rapidity distribution. The third model (Goncalves 
and Machado - GM) describes the photoproduction of 
the vector mesons in UPC events using the QCD color 
dipole approach This model includes nuclear effects 
and parton saturation phenomena. It finds <7 p o=876 mb, 
with a rapidity distribution very different from that of 



the other models. The FSZ and GM models provide pre- 
dictions for the momentum transfer dependence of both 
coherent and incoherent p a production. 

The p° photoproduction on nuclear targets has been 
studied at fixed-target experiments [9j and at RHIC (Rel- 
ativistic Heavy Ion Collider). Previous fixed-target pho- 
toproduction experiments with nuclear targets were done 
at relatively low collision energies @. The PHENIX 
collaboration has studied J/ip photoproduction [l(| in 
heavy ion interactions at RHIC. The Solenoidal Tracker 
at RHIC (STAR) collaboration has published measure- 
ments of the p° production cross-section at a center-of- 
mass energy of v / s/vtv = 130 GeV per nucleon pair [ill ]. 
This work presents results at a higher center-of-mass en- 
ergy of -J snn = 200 GeV. To produce a meson with 
mass m p o at rapidity y p o a minimum photon-nucleon 
center-of-mass energy of W 7 at = [y/SNN m p o exp(y p o)] 1 / 2 
is required. At mid-rapidity this corresponds to about 
12.5 GeV, and in the region \y p o\ < 1 to a range of 
7.6 < W 7 n < 20.6 GeV, which is well above that of 
previous fixed-target photoproduction experiments Q. 

The 7-nucleon collision energy is related to a minimum 
photon energy in the lab frame which is given by 

m„o 

E 1 = -^exp(±y p o) (1) 

The two signs reflect the ambiguity about which nu- 
cleus emits the photon. 

In the rest frame of the target nucleus, the minimum 
photon energy is 27^ times higher than in the lab frame, 
where 7^ is the Lorentz boost of the beam. For ^Jsmn 
= 200 GeV, 7l is about 108; at mid-rapidity, this corre- 
sponds to a photon energy in the target frame of about 
84 GeV. For non-zero rapidities, most of the p° pro- 
duction comes from the solution with the lowest energy, 
so the minimum energy required is even less. To eval- 
uate whether or not the p°s can be produced coher- 
ently, one must compare this energy to the maximum 
photon energy. In the rest frame of the target nucleus, 
this maximum is determined by the uncertainty relation 
E tar 9 et < ^2 _ ^ h / R ^ For 7i = lfJ 8, this maximum 

is approximately 650 GeV. Since the minimal energy re- 
quired to produce a p° is much less than 650 GeV, we 
expect to easily observe coherently produced p° mesons. 

In this study, we have collected data at higher en- 
ergy (s/snn — 200 GeV vs. ^/snn — 130 GeV) and 
with about 10 times more statistics than the previous 
STAR study, allowing for more precise measurements of 
the cross-section. We have extended the previous analy- 
sis by measuring both the coherent and incoherent con- 
tributions to the photoproduction cross-section and we 
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have also measured the spin-matrix elements for p° pro- 
duction. In addition to exclusive p° photoproduction, we 
have studied p° photoproduction accompanied by mutual 
Coulomb excitation, as is shown in Fig. [TJ This pro- 
cess primarily occurs via 3-photon exchange, with one 
photon pro ducing the p°, and one photon exciting each 
nucleus [12L [l3j . Each single-photon reaction is indepen- 
dent, and the cross-sections may be written as an integral 
over the impact parameter 



ct(Au Au -> Au*Au*p° ) = 

= J d 2 b[l - P Bad (b)]P p o(b)P XnA (b)P Xn ^(b) (2) 

where Pu a d{b) is the probability of a hadronic interaction, 
P p o(b) is the probability to produce a p°, and Pxn,i(b) 
and Pxn,2{b) are the probabilities to excite nucleus 1 
and 2 respectively. The three-photon exchange reactions 
are biased toward smaller impact parameters than single- 
photon reactions, leading to a harder photon spectrum 
and an altered rapidity distribution. In mutual Coulomb 
excitation, the nuclei decay primarily by channels involv- 
ing neutron emission. This is attractive experimentally, 
since the neutrons provide simple trigger signals that can 
be detected with the STAR Zero Degree Calorimeters 
(ZDCs) 0. 

One particular nuclear excitation merits special inter- 
est: electroma gne tic excitation to a Giant Dipole Reso- 
nance (GDR) [ig of either one or both ions; GDR in- 
volves particularly low-energy photons. A single GDR is 
the main contribution in the total fragmentation cross- 
section induced by Coulomb excitation in UPCs. GDRs 
usually decay by single neutron emission, which is con- 
sidered to be a major source of beam loss in heavy ion 
colliders Q. 
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FIG. 1: The diagram for p° photoproduction accompanied 
by mutual Coulomb excitation. The latter process proceeds 
by mutual photon exchange; the vertical dashed line shows 
how the photoproduction factorizes from the mutual Coulomb 
excitation. 

The hypothesis of s-channel hclicity conservation for 
vector mesons suggests that the vector meson produced 
in the collision will retain the helicity of the initial pho- 
ton . The differential production cross-section and 
the decay angular distribution of the vector meson can be 
expressed as a function of the vector meson spin density 



matrix, which is represented by the sum of the helicity 
states and encompasses transverse and longitudinal ele- 
ments and their combinations [H. As a consequence of 
parity conservation and the symmetry properties of the 
matrix elements, the decay angular distribution can be 
greatly simplified. If the helicities are conserved in the 
hadronic center-of-mass system, there are only three in- 
dependent helicity amplitudes in the final function [13]. 
In this paper, we report a measurement of these three p° 
spin density matrix elements for px < 150 MeV/c and 
a photon- nucleoli center-of-mass system energy jy » 
10 GeV, beyond previous fixed-target experiments [LSI ]. 



II. EXPERIMENTAL SETUP AND 
TRIGGERING 



This analysis uses data taken with the STAR detec- 
tor at Brookhaven National Laboratory during the 2001 
run. Gold nuclei were collided at ^/snn = 200 GeV 
and the charged particle tracks were reconstructed in a 
cylindrical Time Projection Chamber (TPC) [19]. The 
TPC is a 4.2 m long barrel with a 2 m radius oper- 
ated in a solenoidal magnetic field of 0.5 T. The TPC 
detected charged tracks with pseudorapidity |ry| < 1.2 
and px > 100 MeV/c with an efficiency of about 85 %. 
The TPC is surrounded by 240 Central Trigger Barrel 
(CTB) [20( slats which are plastic scintillator detectors 
spaced every 6 degrees in <£> with complete (hermetic) 
coverage over the full range of rapidity covered by the 
TPC. Two Zero Degree Calorimeters (ZDCs) [3 are sit- 
uated along the beam pipe at ± 18 m from the interac- 
tion point. They have an acceptance close to unity for 
the neutrons originating from nuclear break-up. 

This analysis used data from two triggers: a topology 
trigger and a minimum bias trigger. The topology trigger 
divided the CTB detector into four azimuthal quadrants. 
A coincidence between the left and right side quadrants 
was required, while at the same time the top and bottom 
quadrants were required to be empty. The veto from the 
top and bottom quadrants was used to reduce the trigger 
rate due to cosmic rays. 

The minimum bias trigger required a coincidence be- 
tween the two ZDCs and thus was sensitive to photo- 
production accompanied by mutual Coulomb excitation. 
By eliminating cosmic rays and other extraneous interac- 
tions, this trigger had considerably better selectivity than 
the topology trigger. The ZDCs have sufficient energy 
resolution to count the number of neutrons emitted by 
the outgoing gold nuclei. We distinguish between several 
different excitation modes: XnXn - at least one neutron 
in each of the ZDC detectors, lnln - exactly one neutron 
in each of the ZDC detectors, OnXn - at least one neutron 
in one of the ZDC detectors and none in the other, and 
OnOn - no neutrons in either ZDC. The last two modes are 
only accessible with the topology trigger. A typical ZDC 
spectrum is shown in Fig. [5] In the West ZDC, the ratio 
of ln:2n:3n is 1: 0.48 ± 0.07: 0.42 ± 0.04, while in the 
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East ZDC, we find ln:2n:3n is 1: 0.46 ± 0.08: 0.42 ± 0.04. 
This spectrum allows us to measure the cross-section for 
different excitation states. 
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FIG. 2: ZDC spectra obtained with the minimum bias sample 
after the p° selection cuts are applied, and fit with three Gaus- 
sians. The east ZDC is shown on the left and the west ZDC is 
shown on the right. The ratio of numbers of candidates in the 
West ZDC of ln:2n:3n is 1: 0.48 ± 0.03: 0.42 ± 0.03, while in 
the East ZDC, we find ln:2n:3n is 1: 0.46 ± 0.03: 0.42 ± 0.03. 



III. p° PHOTOPRODUCTION 
A. Event Selection 

This analysis selected events with two oppositely 
charged tracks forming a primary vertex (beam inter- 
action point) and having less than six reconstructed 
charged tracks per event. A p° photoproduction event 
should have exactly two tracks in the TPC, but addi- 
tional tracks may come from overlapping interactions, 
including beam-gas events. The STAR TPC has a 36 ps 
drift time, so any charged particles traversing the TPC 
within ±36 ps may deposit energy which overlaps with 
the tracks of interest. We accounted for the effect of 
these tracks in our analysis by allowing for varying num- 
bers of total tracks in the event, including primary and 
secondary tracks. We analyze only events with exactly 2 
tracks forming a primary vertex. If the cut on the total 
number of tracks is relaxed from 2 to 5, the number of 
reconstructed p° increases by 27 ± 1 %. The results were 
corrected by this factor. 

The reconstruction software formed a vertex from the 
charged particle trajectories. An iterative procedure was 
used to successively remove tracks that were inconsistent 
with the vertex position; after the least consistent track 
was removed, the a new vertex was found. This process 
continued until a vertex was found with an acceptable 
probability. The single track reconstruction efficiency for 
< 1.2 is about 85 %, and the vertex finding efficiency 
for a two-track vertex is 80 ± 2 %. There are several types 
of backgrounds: peripheral hadronic interactions, other 
photonuclear interactions, e + e _ pairs from two-photon 
interactions, and unrelated processes such as beam-gas 
interactions, cosmic ray muons and pile-up events. These 
backgrounds can be reduced by cuts on the total multi- 
plicity, vertex position, and other event characteristics. 



The multiplicity cut suppresses the contribution from 
hadronic and pile-up events. After the cuts on the mul- 
tiplicty, the minimum bias and topology samples con- 
tain 48,670 and 98,112 events, respectively. In order to 
reduce the backgrounds originating from processes like 
beam-gas events, upstream interactions, cosmic rays and 
pile-up events, we selected events with primary vertices 
within 15 cm radially and 100 cm longitudinally (along 
the beam direction) of the center of the interaction re- 
gion. Those two cuts reject approximately 25 % of the 
events. We also required that tracks have at least 14 
hits in the TPC (out of a maximum of 45 possible hits). 
This cut rejected another 30 % of the events. In order to 
retain as many of the incoherently produced p° mesons 
as possible while removing combinatorial background, a 
relatively soft cut on the p° transverse momentum (pt < 
550 MeV/c) was applied. After cuts, the minimum bias 
sample contains 5,011 events, while the topology sample 
contains 14,693 events. 

Backgrounds from two-photon interactions and non-p° 
photonuclear interactions are small. The cross-sections 
for two-photon production of e + e~ in the STAR accep- 
tance are small [2l| . They were a small correction in the 
130 GeV analysis, but the current study requires M„+ n - 
> 500 MeV/c 2 . With this cut, the corrections are negli- 
gible. 

There is also a small background from coherent uj pho- 
toproduction; the ■k + tt~ final state has a negligible 2.2 % 
branching ratio [23], but the 7r+7r~7r° final state has an 
89 % branching ratio. The {it + it~) invariant mass dis- 
tribution from 3-pion decay peaks at a lower value, and 
higher pt than the p° — > 7r+7r~. This is about a 2.7 % 
correction to the incoherent p° cross-section; we neglect 
this here. 

The hadronic interactions produce much higher mul- 
tiplicity final states than the photoproduced p° and can 
be easily distinguished by their total multiplicity [22] . 

Even with the veto from the top and bottom CTB 
quadrants in the trigger, some cosmic-rays remain in the 
topology sample. Particles that pass near the interaction 
region may be reconstructed as a pair of back-to-back 
tracks with net charge 0, net px ~ and y p o w 0. These 
events are removed by applying a cut on the rapidity so 
that the accepted events have |y p o| > 0.01. On the other 
hand, the ZDC energy coincidence requirements largely 
eliminate cosmic-ray contamination in the minimum bias 
sample. 

We use two approaches to estimate the remaining 
backgrounds. As with the 130 GeV analysis, like-sign 
pairs (7r + 7r + and 7r~7r~) provide a good background 
model 11]. That analysis only considered coherent p° 
production; the like-sign background was scaled up by a 
factor of 2.2 to match the data at high px- By definition, 
this treats incoherent p production as a background. 
We use this approach to measure the ratio of directly 
produced 7r + 7r _ pairs to p° production (|-B/A| ratio) for 
the coherently produced p° mesons, since it correctly es- 
timates the combinatorial background. 
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For the rest of the measurements, we use the unsealed 
background in order to retain the incoherent p° signal. 
For incoherent p° photoproduction, we split the invari- 
ant mass histogram into different pt ranges, and fit each 
Pt bin separately to determine the yield. In our fits to 
the M„+„- spectrum the background is parameterized by 
a polynomial. The polynomial function is fixed with pa- 
rameters obtained from the fit of the polynomial function 
to the non-scaled like-sign distribution. These different 
approaches for the background description cause a 3 % 
systematic error in the cross-section measurement. 



B. Efficiency and Acceptance Determination 

The acceptance of the detector was studied using a 
Monte Carlo event generator which is based on the KN 
model d, [T3| to generate events which reproduce the 
kinematic properties and spatial distributions of the p° 
mesons produced via coherent photoproduction. These 
events were passed through a realistic detector simula- 
tion which reproduces detector resolution and efficiency. 
The efficiency includes the detector acceptance, track and 
vertex reconstruction efficiencies, and selection cuts. 

The p° reconstruction efficiency was studied as a func- 
tion of pt, p\i (azimuthal angle in the center-of-mass 
system of AuAu), 8 (polar angle in the center-of-mass 
system of AuAu), y„a (rapidity) and M^-. The mean 
efficiency for minimum bias p°s with \y p o\ < 1 is 44 ± 
2 %. This efficiency is relatively constant with respect to 
Pt and azimuthal angle, but drops as \y p o\ increases, due 
to the TPC acceptance dropping at higher rapidity. The 
mean efficiency for topology-triggered p°s with \y p o \ < 1 
is 11 ± I %, the efficiency drops slowly as pt or \y p o\ 
rises. There is also an azimuthal dependence due to the 
topology trigger veto regions. 

The estimated resolution for pr, y p o and M 7r + 7r - are 
approximately 6 MeV/c, 0.01 and 6 MeV/c 2 respectively 
for track pairs that passed through the p° selection cuts. 



C. Luminosity 

The luminosity for the minimum bias data sample is 
calculated by assuming that the main contribution to the 
total cross-section arises from hadronic production, with 
a known cross-section. The luminosity was measured by 
counting events with at least 14 primary tracks with pt > 
0.1 GeV/c and \y p o\ < 0.5. These events correspond to 
80 % of the total hadronic production cross-section of 
7.2 b [24]. An extra correction is required to remove 
the effects of an unstable dead time caused by the SVT 
(Silicon Vertex Detector). The integrated luminosity of 
the minimum bias sample is measured to be L = 461 
mb _1 with a systematic uncertainty of 10 % which is 
largely due to the uncertainty of the gold-gold hadronic 
cross-section. 



D. Invariant Mass Fit Function 

The invariant mass distribution of track pairs was 
found by assuming that all reconstructed particles were 
pions; no particle identification was needed due to the 
low background level after selection cuts were applied to 
the tracks. The invariant mass distributions for the min- 
imum bias and topology samples are shown in Fig. [3j 

Pion pairs may be photoproduced via an intermediate 
p° , or the photon may fluctuate directly to a ir + ir~ pair. 
The direct process produces a flat M^+k- mass distri- 
bution. The two experimentally indistinguishable pro- 
cesses interfere and the interference is constructive for 



Mir+n- < M„o and destructive for 
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FIG. 3: Top: The invariant mass distribution for the coher- 
ently produced p° candidates from the minimum bias sam- 
ple with the cut on the p° transverse momentum pr < 
150 MeV/c. Bottom: The invariant mass distribution for the 
coherently produced p° candidates obtained from the topol- 
ogy sample with the cut on the p° transverse momentum pr < 
150 MeV/c. The hatched area is the contribution from the 
combinatorial background. The solid line corresponds to Eq. 
[3] which encompasses the Breit-Wigner (dashed), the mass 
independent contribution from the direct tx + -k~ production 
(dash-dotted), and the interference term(dotted). 
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The invariant mass distribution of the p° candidates 
was fit with a relativistic Breit-Wigner [26| function plus 
a contribution for the direct tt + tt~ production and an 
interference (Soding) term [2?], [28|. The background is 
described by a 2nd order polynomial. An estimate of the 
background from the like-sign pairs was used to obtain 
the parameters for the polynomial function. 

The fit function is: 



dN 



where 



A- 



D 



+ u 



(3) 



0.5 (GeV/c) 2 . At mid-rapidity in a collider environment, 
the longitudinal momentum transfer squared from the 
target nucleus t\\ = (m^/2E 7 ) 2 is small (w 2 (MeV/c) 2 ). 
Therefore t « t± ~ p\ and so we can extrapolate the 
ZEUS measurement of the t dependence of l-B/^-l down 
to our average value of t « p T < 0.015 (GeV/c) 2 . Ex- 
trapolating the ZEUS results, we find 0.8, which 
is consistent with our results. The decrease of |-B/A| with 
increasing \t\ and the independence of the polar and az- 
imuthal angles is expected assuming that the form factor 
of the vector meson depends only on t and so no ^/snn 
dependence is expected [25$ . 

In the rest of this paper, we will quote our results in 
terms of momentum transfer squared t. 



r p « = T • (M p o/M T+T -) 



4m 2 )/(M 2 - 4ml)} 



3/2 



(4) 



is the momentum-dependent width of the p° and M p o 
is the mass of the p°, A is the amplitude for the Breit- 
Wigner function, B is the amplitude for the direct tt + h~ 
production and f p is the fixed second order polynomial 
which is used to describe the background. For the mini- 



mum bias data set, T p0 = 0.162 ± 0.007 GeV/c 2 



and M p o 



= 0.775 ± 0.003 GeV/c 2 from the fit. These values are in 
good agreement with the PDG [23| values. The difference 
between the yield obtained with fixed p° width and mass 
position from that obtained without fixing the width and 
mass is about 2 %. Fixing the width leads to an increase 
in the x 2 /DOF of 5 %. Using the above mentioned fit 
procedure, the minimum bias sample contains 3,075 ± 
128 p° candidates, while the topology sample contains 
13,054 ± 124 p° candidates. 

For the minimum bias data, the measured value of 
\B/A\ is 0.89 ± 0.08 (stat.) ± 0.09 (syst.) (GeV/c)" 1/2 ; 
the systematic error due to the background description 
is 3 %. Figure 2] shows that |-B/j4| does not vary signifi- 
cantly as a function of rapidity. Since rapidity is related 
to photon energy (Eq. [T]) this also shows that there is 
no significant variation with photon energy. The same 
energy-independence was seen by the ZEUS experiment 
at HERA |26[ , and is expected in Pomeron exchange. In 
contrast, at lower energies j^j, the direct tt + tt~ compo- 
nent decreased as W rose. The difference may be due 
to p production by meson exchange; the meson exchange 
component is expected to be negligible at RHIC energies. 
The observed ratio is independent of the polar and 

azimuthal angle, as expected. The slight asymmetry in 
the distribution is believed due to differences between the 
two parts of the STAR TPC, and is included within the 
systematic uncertainties. 

Our measured value for is in agreement with 

the 130 GeV STAR result, \B/A\ = 0.81 ± 0.08 (stat.) 
± 0.20 (syst.) (GeV/c)" 1/2 The ZEUS studies 

of 7p -» p°p find \B/A\ = 0.67 ± 0.03 (GeV/c)~ 1 /2 
[26| . but for a slightly different kinematic interval. The 
ZEUS results are for a momentum transfer squared of t < 
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FIG. 4: The ratio |B/j4| as a function of y p o for the minimum 
bias data, obtained by fitting Eq[3]to the invariant mass dis- 
tributions in bins of y p o. 



E. 



Cross-Sections da/dy and d 2 a/dydt for Minimum 
Bias events 



The differential cross-section as a function of rapidity 
for p° photoproduction (da/dy) obtained with the min- 
imum bias sample is shown in Fig [5] The distribution 
for each rapidity bin was fit with Eq. [3] and the p° yield 
extracted. Also shown is a prediction of the KN model; 
the other two models do not include nuclear excitation. 

Fig. [6]shows the p° spectrum as a function of t for min- 
imum bias data. The efficiency correction and luminosity 
normalization have been applied. 

We do not observe the dip in the range 0.01 < t < 
0.015 (GeV/c) 2 predicted by FSZ 0. The p T of the p° 
meson is the vector sum of the photon px and the px 
transferred by the target nucleus. Since the direction 
of the photon and scattering p? are uncorrelated, this 
addition smears out the diffractive dips (29[. 
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The d 2 a /dydt distribution (averaged over \y p o\ < 1) 
is fit to a sum of two exponentials, which correspond to 
coherent and incoherent production: 



d 2 a 
dydt 



= A coh exp (~B coh t) + A inc exp (-B mc t). (5) 



The simple fit function shown in Eq. [5] has two draw- 
backs. The interference between p° photoproduction on 
the two nuclei reduces d 2 a /dydt at small t [H, [29], and, 
in fact, alters the minimum bias t— spectrum at the 20 % 
level for t < 0.01 (GeV/c) 2 . 

A fit over all t values yields exponentials (shown in 
the left column of Table I) that are integrable to give 
the total coherent cross-section (given are the parame- 
ters of the exponentials, not the cross-section). Because 
of the interference at small t, this fit has a poor \ 2 /DOF 
of 79.12/10. A second fit (shown in the right column 
of Table I) is performed over the t range from 0.002 to 
0.3 (GeV/c) 2 . It avoids the region where the interference 
is large and has a y^/DOF of 8.1/7. This fit yields a 
nuclear slope with accuracy comparable to other exper- 
iments. Both fits give similar results for the incoherent 
production. 

The incoherent slope, B inc = 8.8 ± 1.0 (GeV/cT 2 
has not previously been determined in heavy ion col- 
lisions. However, it is comparable to the slope ob- 
served by STAR in dAu collisions [30], and compara- 
ble to the ZEUS result B p = 10.9 ± 0.3 (stat.) ±J;g 
(syst.) (GeV/c)- 2 [26] and HI result B p = 10.9 ± 2.4 
(stat.) ± 1.1 (syst.) (GeV/c)" 2 [p for p° photopro- 
duction on proton targets at comparable t values. The 
HERA data are at higher W 7 n, but the energy difference 
is not expected to introduce a large shift. 

With the second fit region we find the coherent pro- 
duction slope of B coh = 388 ± 24 (GeV/cT 2 , ob- 
tained with the double exponential fit function. For di- 
rect comparison with the previous STAR result, we also 
performed a single exponential fit, which gave -B con = 
363 ± 21 (GeV/c) -2 in agreement with the value ob- 
served at 130 GeV, 358 ± 31 (GeV/c)" 2 [HI- These 
numbers are not directly comparable with fixed-target 
photoproduction data because the photon flux in UPC 
photoproduction is higher on the side of the target near- 
est the photon emitter and lower on the far side of the 
target. The photon flux falls as 1/r 2 , which leads to a 
slightly smaller apparent source size. 

Despite these difficulties, the two exponentials in Eq. 
[5] were integrated analytically to find the total coherent 
and incoherent cross-sections. This approach neglects 
the corrections due to the loss of incoherent cross-section 
when the coherent cross-section is large [32| , but is useful 
for phenomenological comparisons. For \y p o \ < 1, we find 
the ratio 



^incoherent/ ®< 



coherent 



= 0.29 ± 0.03 (stat.) ± 0.08 (syst.) 



for events with mutual excitation (XnXn). 

We have also studied the cross-sections for p° pro- 
duction accompanied by single neutron emission (lnln) 



which is largely due to mutual excitation of Giant 
Dipole Resonances. We did this by fitting the ZDC 
spectra in Fig. [2] and extracting the single neutron 
component from the fit. For \y p o\ < 1, we find 

^"coherent/^helent = 0.18 ± 0.08 (stat.) ± 0.05 (syst). 

The higher crmcoherent/ccohcrcnt for the XnXn sample 
may signal a breakdown of the factorization implicit in 
Eq. [2] because the incoherent p° production transfers 
enough energy to disassociate the target nucleus. This 
effect would lead to additional multiple neutron emission 

M. 




FIG. 5: Coherent p° production cross-section for the mini- 
mum bias data set as a function of y p o (black dots) overlaid 
by the da/dy distribution predicted by the KN model [22l ] 
(solid line). 




0.25 0.3 

t, (GeV/c) 2 



FIG. 6: p production cross-section as a function of the mo- 
mentum transfer squared t, together with the fit of Eq. [5] 
The fit parameters are shown in Table [T] 
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TABLE I: Parameters for the fit to the d 2 o/dydt, Eq. [5] 



Parameter 


t range (0.,0.3) 
(GeV/c) 2 


t range (0.002,0.3) 
(GeV/c) 2 




A coh , mb/(GeV/cy 
B coh , (GeV/c)" 2 


1050 ± 57 


2307 ± 258 


256 ± 12 


388 ± 24 


Amc, mb/(GeV/c) 2 
Bine, (GeV/c)" 2 


21.6 ± 2.4 


24.8 ± 2.5 


7.9 ± 0.9 


8.8 ± 1.0 



Total Cross-Sections 



We have compared three theoretical models to our 
measurements U H S] ■ The comparison is shown in Fig. 

m 

The total production cross-section, do tot I dy, is 
obtained by scaling the mutual excitation results 

and 



XnXn) 



with the scaling factors o-(p° 0n0jl )/(r(p 
v(P°OnXn)/v(P°XnXn) as a function of rapidity. The 
scaling is needed because the efficiency of the topol- 
ogy trigger, which enters into the total cross-section, is 
only poorly known. Therefore the p° production cross- 
section for the events with mutual excitation {XnXn) 
measured with the minimum bias sample was extrap- 
olated based on the ratios a(0n0n)/(r(XnXn)— 7.1 ± 
0.3 (stat.) and o{0nXn)/o(XnXn) =3.5 ± 0.2 (stat.) 
which are measured within the topology sample. Due 
to the limited acceptance in rapidity, we cannot distin- 
guish between the different theoretical models based on 
the shape. However, the amplitude can be used to elim- 
inate models which significantly overestimate the total 
production cross-section in the measured rapidity range. 
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FIG. 7: Comparison of theoretical predictions to the mea- 
sured differential cross-section for coherent p° production. 
The statistical errors are shown by the solid vertical line at 
each data point. The sum of the statistical and systematic 
error bars is shown by the grey band. 

The cross-sections for coherent and incoherent produc- 
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tion for |y p o| < 1 accompanied by nuclear excitation 
are a coh (XnXn, \y p0 1 < 1) = 14.5 ± 0.7 (stat.) ± 1.9 
(syst.) mb and Oi nc (XnXn, \y p o \ < 1) = 4.5 ± 0.5 (stat.) 
± 0.6 (syst.) mb. 

Finding the total cross-sections requires an extrapola- 
tion to the region \y p o \ > 1. This extrapolation is neces- 
sarily model dependent. The KN H2 and FSZ Q calcu- 
lations have similar do tot I dy distributions and so a sin- 
gle extrapolation should work well for them. For the KN 
calculation, the extrapolation factor from cr(|y p o| < 1) 
to a tot is 2.2 ± 0.1 for the events with nuclear break- 
up. We assume that this factor is the same for coher- 
ent and incoherent production. The coherent production 
cross-section extrapolated to the full rapidity range is 
a coh {XnXn,Ml-y) = 31.9 ± 1.5 (stat.) ± 4.5(syst.) mb. 
The total cross-section is 

CTcoh + inc{XnXn, full-y) = 41.4 ± 2.9 (stat.) ± 5.8 (syst.) mb 

the XnXn denotes multiple neutron emission due to nu- 
clear dissociation. 

For p° production accompanied by single neutron 
emission, we find, cr co h(lnln, \y p o | < 1) = 1.07 ± 
0.08 (stat.) ± 0.14 (syst.) mb and eri nc (lnln, |y p o | < 
1) = 0.21 ± 0.09 (stat.) ± 0.03 (syst.) mb. 

The extrapolation factor from \y p o\ < 1 to 47r is as- 
sumed to be the same as that for the XnXn dataset (i.e. 
2.2). The total cross-section for single neutron emission 
is 

CTcoh + inc(lnln, full-y) = 2.8 ± 0.3 (stat.) ± 0.4 (syst.) mb 

Based on the ratio o{p° ' 0n0n ) / 'o(p° XnXn ), we find 
0- coh (OnOn, \y p o\ < 1) = 106 ± 5 (stat.) ± 14 (syst.) mb. 

As with the XnXn data, the extrapolation to 47r is 
model dependent. For the KN model, the extrapolation 
factor is 3.7. For the FSZ model, the factor would be 3.5, 
and for the saturation model GM Q, 2.13. The KN and 
FSZ model factors are similar, what is remarkable since 
the predicted cross sections are differ by 60%. Since the 
KN da/dy matches the XnXn data well, we adopt an 
overall extrapolation factor of 3.7 ± 0.1. With that, we 
find cr coh (0n0n, full-y) = 391 ± 18 (stat.) ± 55 (syst.) 
mb and a total cross-section of events with OnOn (coher- 
ent, incoherent) is 

cr C oh+inc(0n0n, full-y) = 508 ± 24 (stat.) ± 71 (syst.) mb 

It is also possible for a single nucleus to be excited, 
OnXn in this language. The OnXn cross-section includes 
the possibility for either of the two nuclei to dissociate. 
We have checked that we get symmetric results for this 
channel when the signals are in the east or west ZDC; 
those events are added and treated equally. 

This yields the total coherent cross-section 
<7 coh (AuAu -> Au ( *W*y) = 530 ± 19 (stat.) 
± 57 (syst.) mb, and total cross-section (coherent, 
incoherent) 

CT C oh+inc(AuAu -► AuWAu^V ) = 

697 ± 25 (stat.) ± 73 (syst.) mb. 
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Table [HI summarizes the measured coherent and inco- 
herent production cross-sections and compares them with 
results obtained at ^snn = 130 GeV ll]. The measured 
12% increase (with large errors) in the coherent photo- 
production cross-section (going from 130 to 200 GeV col- 
lisions) is much less than is predicted by all 3 models @ , 
[22| and @ , which predict cross-section increases of be- 
tween 70% and 80%. Our results at ^/sWn = 200 GeV/c 
support the //"-nucleon cross-sections used in KN (22|. 

Several sources of systematic error have been consid- 
ered in this analysis. The main sources of the systematic 
errors for the cross-section in the rapidity range \y p o \ < 1 
are 10 % for the luminosity measurement, 3 % due to the 
different approaches for the background description and 
7 % for the applied cuts and fit function. The major addi- 
tional systematic error for the total coherent and incoher- 
ent production cross-sections is 6 % for the extrapolation 
to the full rapidity interval. The error is mainly due to 
the difference between extrapolation factor in KN and 
FSZ models. These uncertainties were added in quadra- 
ture to give the systematic error for the total production 
cross-section. 



G. p" Spin Density Matrix 

The angular distribution of the 7r + and n~ in the 
7— nucleon center of mass frame can be used to determine 
the p° spin density matrix elements. This has previously 
been studied in jp collisions at HERA [26|. There, mea- 
surement of the recoiling proton allowed the jp center of 
mass frame to be determined. STAR does not observe 
the recoiling proton, and so we cannot separate the mea- 
sured pt into contributions from the photon and from 
the nucleon. Furthermore, there is a two-fold ambiguity 
about photon direction. Because of these problems, we 
perform our analysis with respect to the z-axis (beam 
direction). Since the laboratory frame is heavily boosted 
with respect to the 72? center-of-mass frame, this is a 
good approximation; in the target frame, the p° direc- 
tion is within 1-2 mrad of the beam axis. &h is defined 
as the polar angle between the beam direction and the di- 
rection of the 7r + in the p° rest frame. With the two-fold 
ambiguity in photon direction, the +z and — z directions 
are equivalent. This does not affect two terms: cos 2 (0/j) 
and sin 2 (0/ l ), since they are symmetric around 7r/2 (i.e. 
around mid-rapidity). However the term sin(20^) is not 
symmetric around 7r/2 and therefore we are not sensitive 
to the interference between helicity states non-flip to sin- 
gle flip. The azimuthal angle <5>/j is the angle between the 
decay plane and the p° production plane. The produc- 
tion plane of the p° contains the p° and a virtual photon. 
The dependence of the cross-section on &h and cos(0h) 
can be written as follows 1171: 



1 *° = ±.fh 1 _ r 0* ) 

adcos(Q h )d<S> h 4tt 2 V 00 ' 

+ i(3r ° 4 -l)cos 2 (e,) 

-VMe[r^] sin(26 h ) cos($ h ) 
-r 1 i lS in 2 (6 /l )co S (2$ ?i )]. (6) 

The three independent spin density matrix elements 
r 00' r i0' r i-i can be extracted by fitting the two dimen- 
sional angular correlation. The superscripts indicate con- 
tributions from the photon polarization states [3~Ij . The 
element Tqq represents the probability that the p° is pro- 
duced with helicity from a photon with helicity ± 1. 
The element ri 4 _ 1 is related to the size of the interference 
between the helicity non-flip and double flip and SRe^^] 
is related to the interference of non-flip to single flip, 
where Kefr^] stands for the real part of r^. If helicity 
conservation holds, then all three matrix elements will be 
close to zero. 

Figure shows the $h vs cos(0h) correlation fit for 
the minimum bias data set. The measured spin density 
matrix elements are shown in Table IIII1 The method 
used is to fit the invariant mass distributions in bins of 
<f>h and cos(0/i) to determine the yield in each bin. 

The background is accounted for in the fitting function 
with non-scaled like-sign pairs as described in section 
IIII Al The main contribution to the systematic uncer- 
tainty comes from the background subtraction. It was 
estimated by using an alternative approach where the 
scaled invariant mass distribution of the like-sign pairs 
was subtracted from that of the opposite-sign pairs. An 
additional source of systematic error is the uncertainty 
due to the acceptance correction determined by a p° 
Monte Carlo simulation. We estimate the systematic 
error by varying the bin size of <&h and cos(G^). The 
systematic error for the spin density matrix elements is 
obtained by adding the individual uncorrelated contri- 
butions in quadrature. The measured p° helicity matrix 
elements indicate that helicity is conserved within errors 
as expected based on s-channel helicity conservation. 

IV. CONCLUSION 

Photoproduction of p° mesons has been measured in 
the STAR detector at RHIC in Au-Au collisions at ^/snn 
= 200 GeV. Coherent and incoherent p° photoproduction 
has been observed and photoproduction of the p° mesons 
is observed with and without accompanying Coulomb nu- 
clear excitations. The measured increase with energy in 
the total cross-section for photoproduction is much slower 
than proposed in_0] and 0]. However, the Klein and 
Nystrand model j5| is able to describe the data for two 
energy points ^/snn = 130 and 200 GeV. 

The differential cross-section for photoproduction has 
been studied as a function of t, y„o and M^+ n -. The 



11 



TABLE II: The total cross-section extrapolated to the full rapidity range for coherent p production at \/snn = 200 GeV 
accompanied by nuclear breakup and without breakup compared with previous measurements at ^/snn = 130 GeV [111]. The 
first error is statistical, the second is systematic. 
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yJs NN = 130 GeV 


Vsjviv = 200 GeV 


Vsivjv = 200 GeV 


"Lu (mb) 


coherent 


coherent 


coherent + incoherent 


28.3 ± 2.0 ± 6.3 


31.9 ± 1.5 ± 4.5 


41.4 ± 2.9 ± 5.8 


a 0nXn ( mb ) 


95 ± 60 ± 25 


105 ± 5 ± 15 


145 ± 7 ± 20 


^l^in ( mb ) 


2.8 ± 0.5 ± 0.7 


2.4 ± 0.2 ± 0.4 


2.8 ± 0.3 ± 0.4 


vLon ( mb ) 


370 ± 170 ± 80 


391 ± 18 ± 55 


508 ± 24 ± 71 


"total ( mb ) 


460 ± 220 ±110 


530 ± 19 ± 57 


697 ± 25 ± 73 
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FIG. 8: Projections of the two dimensional efficiency cor- 
rected $h vs cos(Oh) distributions obtained with the mini- 
mum bias data set. The solid line shows the result of the 
two-dimensional fit to the data with Eq. [6] and the coeffi- 
cients given in Tab. Mil 



TABLE III: Measured spin density matrix elements compared 
with ZEUS(7p) results. The first error is statistical, the sec- 
ond is systematic. 



Parameter 


Fit result 


7P experiment [26] 


X'fndf 


26/21 




ro°o 4 


-0.03 ± 0.03 ± 0.06 


0.01 ± 0.01 ± 0.02 


«e[rSS] 




0.01 ± 0.01 ± 0.01 


r?i a 


-0.01 ± 0.03 ± 0.05 


-0.01 ± 0.01 ± 0.01 



d 2 a/dydt distribution was fit with a double exponential 



function to isolate the incoherent part of the p° produc- 
tion cross-section. 

The ratio of tt + it~ to direct p° production (|£?/v4|) has 
been studied with respect to polar angle, azimuthal angle 
and y p o; no dependence has been observed as predicted 
in reference [251 ] . 

The Tqq and r'(_ 1 spin density matrix elements for 
the p meson were measured. The small values of Tqq 
and r® 4 L 1 indicate that helicity is conserved within er- 
rors as expected based on s-channel helicity conservation 
(SCHC), and therefore we see no evidence for p° photo- 
production involving spin flip. 
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